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Abstract 

The ability to detect extremely small forces and nanoscale displacements is vital for disciplines 
such as precision spin-resonance imaging [1], microscopy [2j, and tests of fundamental physical 
phenomena [SHS]. Current force-detection sensitivity limits have surpassed 1 aN/\/Hz [6l [7] (atto 
= 10~^®) through coupling of nanomechanical resonators to a variety of physical readout sys- 
tems [U [TUTOj. Here we demonstrate that crystals of trapped atomic ions [IH |T2] behave as 
nanoscale mechanical oscillators and may form the core of exquisitely sensitive force and displace- 
ment detectors. We report the detection of forces with a sensitivity 390ibl50 yN/\/Hz (more than 
three orders of magnitude better than existing reports using nanofabricated devices [7|), and dis- 
criminate ion displacements ^^18 nm. Our technique is based on the excitation of tunable normal 
motional modes in an ion trap [13j and detection via phase-coherent Doppler velocimetry |141 115]. 
and should ultimately permit force detection with sensitivity better than 1 yN/-v/Hz |16| . Trapped- 
ion-based sensors could permit scientists to explore new regimes in materials science where aug- 
mented force, field, and displacement sensitivity may be traded against reduced spatial resolution. 
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Trapped atomic ions exhibit well characterized and broadly tunable (kHz to MHz) normal 
motional modes in their confining potential [ini E] . The presence of these modes, the light 
mass of atomic ions, and the strong coupling of charged particles to external fields makes 
trapped ions excellent detectors of small forces with tunable spectral response [13] . Another 
advantage is that readout is achieved through resonant-fiuorescence detection using only a 
single laser. Previous studies have suggested that by using ions it is possible to measure forces 
approaching the yoctonewton scale, for instance, through experiments on motional heating in 
Paul traps due to fiuctuating electric fields [JBHSO] , or resonant excitation techniques [TTl EI] ■ 
In particular, small forces applied to ions in weak trapping potentials (trapping frequencies 
~0.1 MHz or lower) can excite micron-scale motional excursions resolvable using real-space 
imaging [2HI22]. 

While the intrinsic sensitivity of trapped ions to external forces and fields is well sup- 
ported, it remains an experimental challenge to determine the maximum achievable sensi- 
tivity to a given external excitation as set by systematic limitations including the efficiency 
of a measurement procedure. Establishing ions as components in ultrasensitive detectors 
requires two primary issues to be addressed: a known excitation must be applied to al- 
low precise calibration of the system's response; and it must be possible to compare the 
results of these experiments with the existing literature on detectors based on integrated 
nanostructures. Our aims are to unify the seemingly disparate fields of nanotechnology and 
atomic devices, through use of comparable experimental conditions and a demonstration of 
the potential utility of ion-based sensors in nanomaterials characterization. 

Our system is a crystal of ^Be"^ ions in a Penning trap [TTl [T2l 123] . Ions scatter ultra- 
violet laser light near 313 nm, which is nearly resonant with an internal atomic transition, 
cooling the ions and providing a direct imaging mechanism (see Supplementary Information 
for details). Scattered fiuorescence is detected using an imaging system connected to a CCD 
or photomultiplier tube. For these experiments we focus on ion crystals [23H2B] with n~100 
ions in a two-dimensional planar array (perpendicular to the direction of the detection laser), 
having a diameter ~300 /im. 

In our detection technique, known as laser Doppler velocimetry [IH[l5|, uniform ion mo- 
tion at the axial center-of-mass (COM) frequency (set at u;^/27r = 867 kHz for this work) 
parallel to the propagation direction of the cooling laser beam modulates the intensity of 
resonant fiuorescence due to Doppler shifts. Under oscillatory ion motion, and with the laser 
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detuned from resonance near the Doppler-cooling-efficiency maximum [15] the intensity of 
ion fluorescence is modulated at uz with amplitude approximately proportional to the mo- 
tional amplitude (Fig. la). For uniform light collection from all ions and detection laser 
power below saturation, the total detection rate of scattered photons is np [1 + {2/'^)kzcoM], 
where n is the ion number, p captures all hardware parameters including quantum efficiency 
of the detector, the intensity of illumination, etc., k is the wavevector of the Doppler detec- 
tion laser parallel to the ion motion, and zcom is the velocity of the axial center-of-mass 
coordinate. 

We consider an impulse-style measurement in which a force, sin(ci;rft) is applied to the 
ions for a fixed drive time, td ^ l-n jbJd^ during which the Doppler detection laser is turned 
off. As such, radiation damping during the excitation does not occur. We define to be 
the amplitude of the total force applied to a crystal of n ions, and i^J*""-* to be the amplitude 
of the force applied to a single ion of mass m. After application of this drive pulse, an ion 
crystal of n ions will undergo a steady-state sinusoidal oscillation with velocity 



and the oscillation phase is 



assuming no damping in a harmonic confining potential, mcj^-ZcoM/^- Near resonance 
{ujz ~ Wd), Eq. 2 can be written v = {Fd/2nm)td. 

Ion motion induced by the application of this driving force will be superimposed on a 
noisy background due to the finite temperature of the axial mode and stray electric fields. 
The ability to discern a small oscillating signal of known frequency in the presence of a 
large background is well established through use of phase-sensitive detection [27]. In these 
techniques, shown to be applicable in almost any setting, a signal of interest is discriminated 
through synchronization of a target system's response to a master oscillator that produces 
or modulates the excitation. In doing so, broadband environmental noise is excluded, and 
only the integrated noise over the narrow bandwidth of the measurement is germane. 




(1) 



where for jcj^ — w^l /w^ ^ 1, the amplitude v is given by 




(2) 
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The temporal modulation of ion fluorescence due to the applied drive is detected phase- 
coherently by recording scattered-photon arrival times using a photomultiplier tube relative 
to a trigger synchronized to the external drive force. Photon arrival times relative to the 
drive force are then determined using a time-to-amplitude converter (TAG) or multichannel 
scaler over iterations of the experiment. Any noise, however, is not phased with the drive, 
and averaging A^-times increases the contrast of the desired signal relative to background 
noise by a factor of -y/iV. 

Precise characterization of the the achievable force- and displacement-detection sensitiv- 
ity requires the application of a well calibrated drive force. Such a force is generated using 
the electric fleld from an RF voltage applied to an endcap electrode on the trap, similar to 
resonant electric- fleld excitation techniques used in cantilever experiments [28] (Fig. lb, c). 
The electric fleld at the location of the ions is calibrated using a measurement of the static 
deflection of a planar ion crystal by sideview imaging under application of a static voltage to 
the same endcap electrode (see Supplementary Information) . As an example, for an applied 
zero-to-peak voltage of 165±10 /iV (nominal RF power, Prj? = -70 dBm), we determine an 
electric fleld of 1.8±0.1 mV/m at the location of the ions, and a corresponding force Fq"""^ = 
290±18 yN per ion. 

Figure Id shows a typical histogram of the arrival times of the flrst detected photon syn- 
chronized to the drive for = ooz- For short times no photons are detected due to hardware 
delays (mainly the response of the acousto-optic modulator (AOM) switch). Once photon- 
detection events begin accumulating after ~4 /is, we flnd a bunching of photon arrival times 
with a period commensurate with the 1.15 fis period of the COM oscillation. Phase infor- 
mation is captured in the absolute locations of histogram maxima along the time axis. This 
approach allows the target system to have a variable shot-to-shot phase difference from the 
master oscillator as experimental parameters are swept, thus providing more information 
than standard lock-in detection in which the detection phase offset is flxed. 

We study the Doppler velocimetry signal as a function of Ud and td- The measurement 
scans the drive frequency and records a histogram of stop-pulse delays relative to the drive- 
force trigger (the start pulse). The modulation of the scattered photon rate is plotted as a 
colorscale, after correcting for an exponential decay factor due to the triggering technique 
(see Supplementary Information). The flrst column of Figure 2 (Fig. 2a, d, g, j) shows 
experimental measurements of COM-mode excitation, while the second column (Fig. 2b, e. 
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h, k) shows theoretical calculations based on the formulae above. In these calculations we 
fix td and Uz and allow a variable drive strength and a background, both of which are held 
the same for all t^. 

Theoretical calculations match well with experimental data, replicating both qualitative 
and quantitative features. Doppler velocimetry indicates a linear phase shift in the oscil- 
lator as Ud is tuned through resonance, and the excitation of the oscillator is zero when 
~ ^z\ = '^'^ /td, such that during the driving period the external force and system re- 
sponse desynchronize and resynchronize [22] • Moreover, in addition to the central resonance 
feature, oscillation sidelobes appear, separated by the detuning period 2n/td, and we have 
seen ten sidelobes for strong RF excitation. In the experimental data we also observe a 
damping of the oscillation strength as a function of the delay time due to radiation pressure 
from the detection laser. This effect is not accounted for in the calculations shown in the 
second column of Fig. 2. 

The linewidth of the resonance scales as the inverse of td (Fig. 2c, f, i, 1). Here we plot 
the oscillation amplitude (using the standard deviation of the photon-arrival-time histogram 
for all t as a proxy), as a function of Ud, and find that the amplitude of the measured ion 
velocity matches well with theoretical predictions over all Ud- The quality factor of the COM 
mode is limited by the presence of dark ions produced by background gas collisions and by 
power-supply instabilities, but is high relative to the Fourier-limited linewidth of uz for the 
values of td shown here. Agreement between data and theory breaks down for the largest 
values of td which induce the largest ion velocities. Under these conditions, the ions are 
driven out of the linear-response regime of Doppler velocimetry, or even to the blue side of 
the resonance, resulting in reduced scatter rates manifested as dips in the central oscillation 
lobe (Fig. 2j, 1). Reducing the RF excitation removes these effects. The nonlinear velocime- 
try response can also be used as an independent calibration of the applied force (via Eq. 2 
near resonance), agreeing to within a factor of order unity with the electric-field calibration 
described above (see Supplementary Information) . 

We analyze the force detection sensitivity by systematically reducing the excitation 
strength and examining the system response through Fourier analysis. Fig. 3a shows a 
one-dimensional slice of the time-domain Doppler velocimetry signal for diminishing values 
of Fd, accomplished by reducing P^p (td is fixed at 1 ms) using a planar array of n =130±10 
ions. In order to minimize the number of required experimental cycles we employ a mul- 
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tichannel scaler and collect more than one photon (~3-5) per experimental cycle. 
is varied from F^^"^^ to Fo^'°"^/100, and the magnitude of the system response is reduced 
commensurately. We Fourier transform these data (omitting the systematic delay for short 
times) and display the results on a semilog plot in Fig. 3b, demonstrating a spectral peak 
near the COM resonance frequency with diminishing signal-to- noise ratio (SNR). 

With F^'"''^ = Fo^'°"VlOO = 2.9±0.18 yN/ion and n =130±10, the total force on the 
array is — 377±37 yN. From the spectral peak at uz we calculate the SNR (here 2.3) 
and account for the measurement bandwidth to extract a force-detection sensitivity ~ 1,200 
yN/v^Hz for SNR = 1 (see Supplementary Information). 

The COM axial amplitude of an excitation with may be expressed as 

ZCOM = 7^ ■ (4) 

We characterize sensitivity to spatial displacement by considering small values of both Fd 
and td- For an excitation (as above) with td =1 ms and Fd = 377 yN and n =130, we are 
able to discriminate spatial displacements of Zcom ~ 18 nm (sensitivity ~ 58 nm/v^Hz). 
This amplitude is comparable to the thermal axial extent of the COM mode for a crystal 
of 130 ions (11 nm for T — 0.5 mK), and an order of magnitude less than the thermal 
axial extent of an individual ion in the array (120 nm). These absolute and relative sizes 
indicate that real-space imaging of ion excitation is not possible using, e.g. a CCD camera. 
An excitation of this magnitude is also comparable to common values associated with 
nanomechanical resonators. 

Characterization of the minimum achievable force detection sensitivity is accomplished 
by producing a large motional response to a given applied force. For a given f'^"^^ the 
Doppler velocimetry signal near resonance grows linearly with td, but detection sensitivity 
only increases (gets worse) as ^Jtd■ We have confirmed this scaling using ion crystals oin — 
60, 130, and 530 ions, observing an approximately linear increase in measured SNR up to 
td = 10 ms (without bandwidth normalization). For larger values of td the experimental 
drifts of cjz on timescales comparable to the data acquisition time become significant 
relative to the narrowed Fourier response of the driven oscillator. 

The smallest force we detected is ~170 yN using a crystal of n = 60 ± 5 ions and 
td —10 ms, as described above. Averaging over data acquired with different values of Fd 
(normalizing by the SNR for different values of Fd should yield similar values of force 
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detection sensitivity) we find a minimum force detection sensitivity of 390±150 yN/A/^iz. 
Tlie experimental uncertainty includes statistical fluctuations in averaging, uncertainty in 
the ion number, uncertainty in our calibration of the applied electric field, and imprecision 
in the calculation of the SNR. The improvement in this value relative to that extracted 
from Fig. 3b is derived from both the increased drive time and the reduced ion number 
(see Supplementary Information) . 

It is important to distinguish between the measured systemic detection sensitivity 
(force or displacement) and the intrinsic sensitivity of the ions. The ions respond to the 
external stimulus in a bandwidth set by t^^, for these experiments >3000 times wider 
than T]^/. Thus, the intrinsic ion sensitivity is at least 50 times better than reported in 
our measurement, and in line with sensitivities that may be derived from heating rate 
measurements [TSVEU] . However, in any experiment it is necessary to extract information 
from the system and account for inefficiencies in the readout technique when calculating 
the system's sensitivity. 

Considerable systematic improvement is possible before our system becomes limited by 
the intrinsic sensitivity of trapped ions to forces and fields. Realistic experimental im- 
provements may be achieved by increasing the measurement bandwidth, primarily through 
improved light-collection efficiency (measurements are currently shot-noise limited). For 
instance, experimental modifications allowing for large-solid-angle light collection and 
sampling during measurement at the Nyquist limit both reduce the required number of 
experimental cycles, A^, and hence tm- Such improvements could allow realization of force- 
detection sensitivity ~1.7 yN/A/Hz, comparable to previously published calculations [16] 
(see Supplementary Information). Additionally, systematic improvements that increase 
the stability of uz, will permit longer drive times and hence improved force-detection 
sensitivities; our system is currently limited by the effects of background gas collisions and 
power-supply instabilities. 

We have successfully employed this technique for the detection of optical dipole forces 
induced by Raman lasers. Additionally, the Penning trap permits experiments using well 
over n ^ 10^, thus allowing optimization of sensitivity to electric fields (large total charge 
desirable) rather than applied forces (small total mass desirable). Given a crystal of this 
size and similar experimental conditions we could likely achieve an electric field sensitivity 
of ~500 nV m~^/A/Hz, with intrinsic sensitivity limits still lower. 
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The experiments we have presented have demonstrated phase-coherent excitation and 
detection of yN-level forces and nm-scale displacements induced by oscillating electric fields 
using trapped ions. Our measurements have validated published calculations [IB] suggesting 
that force-detection sensitivities of ~1 yN/-\/Hz are possible for single ion experiments, 
owing to the light mass of harmonically bound trapped ions and the presence of a strong 
readout technique. Vitally, we have included in our detection-sensitivity measurements all 
relevant readout times, experimental dead times, hardware delays, and the like, making 
this a very conservative estimate of achievable force-detection sensitivity. 

Beyond simply measuring extremely small forces, the detection technique we employ 
provides an ability to discriminate motional excitations more than an order of magnitude 
smaller than the thermal extent of a single ion, and deep below the resolution-limits imposed 
by typical imaging systems. Because the size of an ion's motional excursion scales as co^^, 
for forces of the magnitude we detect it is difficult to directly image motional excitations 
at frequencies larger than ~0.1 MHz. Thus phase-sensitive detection of ion motion may 
enable high-frequency measurements as a part of a tunable, broadband detector. 

One potential shortcoming of this detection mechanism relates to the use of a phase- 
synchronous protocol. In circumstances such as detector calibration or the application of 
an optical dipole force it is straightforward to synchronize detection events to the externally 
applied force. In general the problem becomes significantly more complicated when trying 
to measure an unknown force. However, it is typically possible to modulate an unknown 
external signal near uz in such a way as to permit synchronization of detection events ■ 
For instance, in surface science or field mapping studies such as proposed in [16], a drive 
signal in the target device (e.g. a current bias) could be modulated directly, producing 
a time varying, detectable force or field. Synchronization with the modulation envelope 
signal would then permit performance similar to that achieved in our experiments. In 
such an approach the use of an ion-based detector could provide a means to map forces 
and fields at target-probe distances not generally accessible using cantilevers or other 
comparable approaches. As an additional example, in studying the dynamics of surface 
charge accumulation, a laser producing photo-excited charge could be modulated in order 
to provide a synchronization signal [30]. Doing so, and then using phase-synchronous 
detection, would permit a study of high-frequency charging dynamics that are not accessible 
in standard "time-integrated" detection techniques. 
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The results we have presented build on a strong foundation laid by previous trapped-ion 
experiments, and suggest that ion-based sensors may form a vital tool for sensitive 
measurements in nanoscale science. Their utility will be maximized in applications where 
lateral spatial resolution may be traded against increased overall sensitivity and the ability 
to detect forces and fields with standoff distances greater than achievable using detectors 
based on integrated nanodevices. Reahstic, field-deployable sensors will likely require the 
use of integrated ion-trap chips with asymmetric potentials for directional force detection 
and three-dimensional spatial mapping. It may also be desirable to exploit sub-Doppler 
cooling and sideband-detection mechanisms [TTf ITS] for the measurement of stochastic fields 
and forces. 
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FIG. 1: Phase-coherent Doppler velocimetry. (a) Schematic of the atomic resonance employed for detection of ion motion, 
7/27r ail9 MHz as a function of the detection laser detuning from the atomic resonance frequency. The detection laser beam 
is oriented perpendicular to the plane of a 2D ion crystal (Inset). An oscillating ion array periodically traces out a path on 
the atomic resonance profile schematically illustrated by the solid line. The extent of the excursion is set by the magnitude 
of maximum ion velocity. Vertical dashed lines indicate the range of Doppler shift associated with a given excursion along 
the resonance profile, (b) Schematic of the Doppler-detection system based on photon-arrival-time measurements. Disc. 
= Discriminator, TAG = Time-to-Amplitude Converter, MCS = Multi Channel Scaler, PMT = photomultiplier tube, (c) 
Schematic of pulse sequencing/triggering for phase-coherent detection. This scheme is based on previous studies of micromotion 
nulling in a Paul trap |14l and studies of plasma oscillations in a Penning trap |15| . with the important distinction that the 
excitation and detection are segregated into different parts of the measurement procedure. = oscillating drive force due to 
the ac electric field, with dotted line showing the pulse envelope defining the pulsed-excitation period. Start pulses synchronous 
with the Drive are output from a discriminator and fed into the TAC or MCS. A Gate pulse ensures that only the last Start 
pulse of the excitation period triggers the TAC/MCS to begin data acquisition. The Detection Laser is turned on using an RF 
switch controlling an Acousto-Optic Modulator, after which scattered photons may be detected. Stop Pulses are generated by 
the detection of scattered photons using the PMT and are fed to the TAC/MCS. (d) Histogram of photon arrival times relative 
to start-pulses generated synchronously with an RF drive of the COM mode on resonance. Horizontal axis corresponds to the 
time delay, t indicated in panel (c). Photon arrivals are bunched with periodicity given by the driven COM oscillation period, 
following hardware delays. Solid blue line is an exponential fit to the data used to remove a background scattering rate (see 
Supplementary Information). 
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FIG. 2: Phase-coherent detection of the COM mode by RF excitation. _p^'°"' = F^^^°"\ 
to exponential decay in photon arrival times for four driving pulse durations, expressed as a colorscale for experiment (a, d, g, 
j) and theory (b, e, h, k). Horizontal axis represents arrival delay from start pulse and vertical axis represents drive frequency. 
Amplitude of oscillations decays with delay time due to radiation damping from the detection laser and is not accounted for in 
theoretical plots. Third column: standard deviation of photon arrival times as a function of drive frequency. Each data point 
represents the standard deviation of a horizontal slice of the two-dimensional plots (left) and illustrates resonant excitation 
of the COM mode as in Fig. Ic (standard deviation over the whole measurement period is used as a proxy for oscillation 
amplitude). Solid lines represent theoretical fits using fixed drive times with force strength and a constant offset used as free 
parameters. Each row of plots in the figure corresponds to a fixed drive time. Fit parameters extracted from 200 drive 
duration used for longer drive periods. Breakdown in fit quality on resonance increases with drive time as strong excitation 
leads the ions to be shifted out of the linear response regime to the blue side of the Doppler resonance, hence decreasing 
fluorescence. 
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FIG. 3: Calibration of force-detection sensitivity by Fourier analysis, (a) Temporal response to applied = nF^*°"' with 

n = 130 for decreasing drive strength and = 1 ms. Total experiment time tm Ri56 s for each trace including all measurement 
and dead time (recooling, photon detection, hardware delays), corresponding to ^40,000 excitation/detection cycles. Observed 
decay of oscillation magnitude due to radiation damping occurring once the detection laser is turned on. (b) FFT of temporal 
response traces recorded in panel (a) above, plotted on a semilog scale, with same color-coding as in panel (a). Spectral peak 
apparent at tocOM diminishes in strength with decreasing p^°^\ Spectral peak has SNR Ki 2.3 for = O.OIOFq*""^ (see 

text). 
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